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Abstract The Aharonov-Bohm (AB) problem for vector 
bosons by the Duffin-Kemmer-Petiau (DKP) formalism is 
analyzed. Depending on the values of the spin projection, 
the relevant eigenvalue equation coming from the DKP for¬ 
malism reveals an equivalence to the spin-1 /2 AB problem. 
By using the self-adjoint extension approach, we examine 
the bound state scenario. The energy spectra are explicitly 
computed as well as their dependencies on the magnetic 
flux parameter and also the conditions for the occurrence 
of bound states. 
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1 Introduction 

The Aharonov-Bohm (AB) effect [1] has been an usual frame¬ 
work for investigating the arising of phases in the wave func¬ 
tion of quantum particles in various physical models and 
has inspired a great deal of investigations in recent years. 
In the AB effect, the vector potential due to a solenoid gains 
an extraordinary physical meaning. It can affect the quan¬ 
tum behavior of a charged particle that never encounters 
an electromagnetic field. This phenomenon is intimately re¬ 
lated to a non-local boundary condition which relates the 
change in the phase of an electron wave function to the 
amount of flux in the solenoid. The interest in this issue ap¬ 
pears in the different contexts, such as solid state physics [2], 
cosmic strings [3-11] K'-Poincare-Hopf algebra [12, 13] , 
5-like singularities [14], supersymmetry [15], condensed 
matter [16], Lorentz symmetry violation [17-19], quantum 
chromodynamics [20], general relativity [21], nanophysics 
[22], quantum ring [23-26] , black hole [27] and noncom- 
mutative theories [28-30]. 
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In the AB problem of spin-1 /2 particles a two-dimensional 
5-function appears as the mathematical description of the 
Zeeman interaction between the spin and the magnetic flux 
tube [14]. This interaction term is known to cause a split¬ 
ting on the energy spectrum of atoms depending on the spin 
state. In AB problem of spin-1 particles [31, 32], however, 
this characteristic is also present. In Ref. [31], where the 
authors address the AB problem for spin-1 Yang-Mills par¬ 
ticles, it was established that, for the case of spin-1 /2, qua¬ 
sibound states exist for all noninteger flux parameter. The 
existence of these states is related to the penetration of the 
magnetic flux tube by the particle, which is sufficient to pro¬ 
duce sensitivity to the sign of the flux. The difference for 
spin-1 Yang-Mills particles is that the quasibound states ex¬ 
ist only for discrete values of the magnetic flux tube, so that 
penetration occurs only for flux values in a set of measure 
zero. 

In this work, we solve the spin-1 AB problem for bound 
states in the context of the DKP formalism. In our approach, 
we consider the idealized picture of a magnetic flux tube of 
null radius which allows the particles to access the r = 0 re¬ 
gion in a controlled way. Unlike the approach taken in Ref. 
[31], here, we modulate the problem with general bound¬ 
ary conditions. When the spin projection = 0, the radial 
operator can be expressed as a modified Bessel differential 
equation. In this case, the system does not admit bound- 
sate solutions. On the other hand, when the spin projection 

= 1,-1, as we mentioned above, we have the presence 
of a 5-function potential in the equation of motion. As is 
well-known in quantum mechanics, the 5-function poten¬ 
tial guarantees at least one bound state for the particle and 
this property is independent of its spin. For the system con¬ 
sidered here, in first sight, the inclusion of the spin projec¬ 
tion element leads to an equation of motion equivalent to 
the equation for the spin-1 /2 AB problem. Because of this, 
the problem can be addressed by the self-adjoint extension 
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method [33, 34] with the application of appropriate bound¬ 
ary conditions. After imposing the boundary conditions, we 
finally determine the bound states of the vector bosons in 
terms of the physics of the problem, in a very consistent way 
and without any arbitrary parameter. 

2 A short review on Duffln-Kemmer-Petiau equation 

The first-order DKP formalism [35-38] describes spin-0 
and spin -1 particles and has been used to analyse relativis¬ 
tic interactions of spin -0 and spin -1 hadrons with nuclei 
as an alternative to their conventional second-order Klein- 
Gordon (KG) and Proca counterparts. Although the formalisms 
are equivalent in the case of minimally coupled vector in¬ 
teractions [39-41], the DKP formalism enjoys a richness 
of couplings not capable of being expressed in the KG and 
Proca theories [42,43] . Indeed, the DKP formalism has been 
widely used in the description of many processes in elemen¬ 
tary particle and nuclear physics and it proved to be better 
than the KG formalism in the analysis of A 73 decays, the 
decay-rate ratio F(t 7 —> 77 )/F( 7 r*’ —> 77 ), and level shifts 
and widths in pionic atoms [AA 4 6]. The DKP formalism 
has also applications in other contexts, as such, in noncom- 
mutative phase space [47], in Very Special Relativity (VSR) 
symmetries [48], in Bose-Einstein condensates [49, 50], in 
topological defects [51], in thermodynamics properties [52], 
in topological semimetals [53], in noninertial effect of rotat¬ 
ing frames [54], among others. 

The DKP equation for a free charged boson is given by 
[38] (with units in which Ti = c = 1) 

= ( 1 ) 

where the matrices satisfy the DKP algebra 

( 2 ) 

and the metric tensor is = diag (1, — 1, — 1, — 1). That al¬ 

gebra generates a set of 126 independent matrices whose ir¬ 
reducible representations are a trivial representation, a five¬ 
dimensional representation describing the spin -0 particles 
and a ten-dimensional representation associated to spin -1 
particles. The DKP spinor has an excess of components and 
the theory has to be supplemented by an equation which al¬ 
lows to eliminate the redundant components. That constraint 
equation is obtained by multiplying the DKP equation by 
1 — namely 

ij3^j3‘’j3“5,'R = M(l(3) 

where j runs from 1 to 3. This constraint equation expresses 
three (four) components of the spinor by the other two (six) 
components and their space derivatives in the scalar (vec¬ 
tor) sector so that the superfluous components disappear and 


there only remain the physical components of the DKP the¬ 
ory. The second-order KG and Proca equations are obtained 
when one selects the spin-0 and spin-1 sectors of the DKP 
theory. A well-known conserved four-current is given by 

yM = (4) 

where the adjoint spinor *¥ is given by "R = with 

77 ** = 2j5^j5^ — 1 in such a way that (the 

matrices are Hermitian with respect to rj^). Despite the 
similarity to the Dirac equation, the DKP equation involves 
singular matrices, the time component of given by (4) 
is not positive definite and the case of massless bosons can 
not be obtained by a limiting process [55]. Nevertheless, the 
matrices plus the unit operator generate a ring consis¬ 
tent with integer-spin algebra and 7® may be interpreted as a 
charge density. The factor 1 /2 multiplying of no im¬ 

portance regarding the conservation law, is in order to hand 
over a charge density conformable to that one used in the 
KG theory and its nonrelativistic limit [56]. 

3 Interactions in the Dnffln-Kemmer-Petiau equation 

With the introduction of interactions, the DKP equation can 
be written as 

{ili^d^-M-U)^ = 0, (5) 

where the more general potential matrix U is written in terms 
of 25 (100) linearly independent matrices pertinent to five 
(ten)-dimensional irreducible representation associated to 
the scalar (vector) sector. In the presence of interaction, 
satisfies the equation 

+ = ( 6 ) 

Thus, if U is Hermitian with respect to then four-current 
will be conserved. The potential matrix U can be written 
in terms of well-defined Lorentz structures. For the spin- 
0 (scalar sector) there are two scalar, being two vector and 
two tensor terms [42], whereas for the spin-1 (vector sector) 
there are two scalar, two vector, a pseudoscalar, two pseu¬ 
dovector and eight tensor terms [43]. The condition (6) has 
been used to point out a misleading treatment in the recent 
literature regarding analytical solutions for nonminimal vec¬ 
tor interactions [57-59]. 

3.1 Duffin-Kemmer-Petiau equation with minimal 
electromagnetic coupling 

Considering only the minimal vector interaction, the DKP 
equation for a charged boson with minimal electromagnetic 
coupling is given by 

(i/3^D^-M)>F = 0, 


(7) 
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where the covariant derivative is given by = 9^ + wA^. 
In this case, the constraint equation (3) becomes 

= M (1 - jSOjgO) (8) 

and the four-current 7^ retains its form as (4). 

3.2 Vector sector 

Now, we discuss the vector sector (spin-1 sector) of the 
DKP theory. To select the physical component of the DKP 
field for the vector sector (spin-1 sector), we define the op¬ 
erator [60] 

= (9) 

which satisfies R^'' = and R^'' = —R''^. Moreover, 

as it is shown in Ref. [60], and transform as a 

(pseudo)vector and (pseudo)tensor quantities under an in¬ 
finitesimal Lorentz transformation, respectively. From the 
above definitions, the following property is obtained; 

R^yp^=Ri^gycc_RV^,ia (jQ) 

In this way, by applying the R^ and R^'’ operators to the 
DKP equation (7), we obtain 



(11a) 

^ ^ M ’ 

(lib) 


(11c) 

which leads to 


+M^ {R'''P)=0, 

(12a) 


(12b) 

where = df^Ay — dyA^. These results tell us that all el- 


ements of the column matrix obey the Proca equation 
interacting minimally with an electromagnetic field. So, this 
procedure selects the vector sector of DKP theory, making 
explicitly clear that it describes a spin-1 particle embedded 
in a electromagnetic field. 

According to Ref. [61], we can rewrite Eq. (12) in the 
form 


Following the same procedure of Ref. [61], Eq. (13) be¬ 
comes 

[Df,D^+M^-e{S-B)]R'F = 0. (14) 

with B = V X A and the spin operator S = (5',5^,5^) is 
expressed by 

/O 0 0 \ / 0 0 /\ /O -i 0\ 

= 0 0 , 5^ = 0 0 0 , 5^ = i 0 0. (15) 

Vo / 0/ X-iQoJ \q 0 0/ 

In this stage, it is worthwhile to mention that the last term in 
equation (14) is the called Pauli term, which is crucial to give 
meaning to the term that explicitly depends of the spin. Also, 
we can mention that the Pauli term is only important for the 
vector sector (spin-1 sector) of the DKP theory, because that 
term is absent for the scalar sector (spin-0 sector). For this 
reason, we only focus the vector sector of the DKP theory. 

If the terms in the potential = (Aq, A) are time-independent 

one can write 

'F(r,f) = v/(r)e-'''V (16) 


where E is the energy of the vector boson, in such a way 
that the time-independent DKP equation for the vector sec¬ 
tor becomes 


{p-eAf+M^ 


{E-eAof 


e(S-B) 


RV/ = 0. 


(17) 


In the next section, we apply Eq. (17) to AB problem, giving 
a focus to spin effects through the e (S • B) term. We shall see 
later that by carefully modulating the radial operator (17) 
with boundary conditions, it can provide both bound and 
scattering states. However, we emphasize only bound states. 


4 The Aharonov-Bohm problem 

Let us consider the particular case where the boson moves 
in the presence of the AB potential (Aq = 0). The vector 
potential in the Coulomb gauge is 

0 

eA = --^, (18) 

P 

where 0 is the flux parameter. The potential in (18) provides 
a magnetic held perpendicular to the plane (p, (p), namely 


-f M“] R'''¥-D''D^R^'¥ -^■^R'' = 0, 

(13) 

where S^'' = The term is called the 

anomalous term because it has no equivalent in the spin- 
1 /2 Dirac theory [38]. However, it has been shown in Refs. 
[39, 40] that such an anomalous term disappears when the 
physical components of the DKP held are selected. 


eB = -0^1^z, (19) 

P 

where B is the magnetic held due to a solenoid. If the solenoid 
is extremely long, the held inside is uniform, and the held 
outside is zero. However, the boson is allowed to access the 
p —0 region. In this region, the magnetic field is non-null. 
If the radius of the solenoid is po ~ 0, then the relevant mag¬ 
netic field is B ^ 5{p) as in (19). 
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4.1 Aharonov-Bohm problem for the spin-1 sector 


Now, we consider the effect of Aharonov-Bohm flux field 
on vector bosons. Substituting Eq. (18) in Eq. (17), we ob¬ 
tain 


0 - 

tV -f --(p 
i p 


-(j)S 


m 


R\j/= -M^)R\I/, 


( 20 ) 


where S is the matrix 



( 21 ) 


and 

„ . d d 

V — P3-h -^*3“’ 

dp p d(p dz 


( 22 ) 


is the gradient operator in cylindrical coordinates. At this 
stage, we can use the invariance under boosts along the z- 
direction and adopt the usual decomposition 

Rxtf = R‘xi/{p,(p,z)=fl;Hp)e‘"‘’^e‘P^^, (23) 


with m € Z. Inserting Eq. (23) into Eq. ( 20), we get 
D/i')(p)=^Vm'(p), (24) 

with k = —M^ — p^, where 

D = Do-f0s'-®, (25) 

P 


1 d (m + <j))^ 
dp^ p dp p^ 


(26) 


and s' = (1, — 1,0) represents the eigenvalues of the operator 
S acting on the DKP spinor fm (p). Equation ( 24) describes 
the quantum dynamics of vector bosons in the presence of 
the Aharonov-Bohm potential. Erom (24) we can see that 
scattering states occur only if A: G R, whereas bound states 
occur only if ^ = i|^|. 

At this level, it is worthwhile to note that the solution for 
this problem can be separated in two cases. The first case is 
when the spin projection = 0. In this case the Pauli term 

C3't 

is absent and the radial operator D becomes Do and fm ^ can 
be expressed as a solution of the modified Bessel differential 
equation. We can see that the system does not admit bound- 
state solutions. On the other hand, the second case is when 
i * = 1,-1. In this case the radial operator D is equivalent 
to Eq. (29) of Ref. [62] (see also Refs. [14, 63, 64]) which 
governs the quantum dynamics of the usual spin-1 /2 AB 
problem, namely 


dP- 


1 d {m + (j))^ 
r dr 


+ 0 S 


8{r) 


(27) 


where s is twice the spin value, with s = +1 for spin “up” 
and s = — 1 for spin “down”, and 0 is the flux parameter. 

In order to study the dynamics of the system for = 
1,-1, it is necessary to solve Eq. (24). However, as for the 
case studied in Ref. [62], it also involves a singularity in 
the r = 0 region. The appropriate method for studying this 
problem is the self-adjoint extension method of operators 
in quantum mechanics. By exploiting the nature of the 5- 
function in Eq. (25), we can see that the system admits at 
least a bound state. In addition, we also must verify that 
the operator Do in Eq. ( 26) admits or not self-adjoint ex¬ 
tensions. It is known that when we consider Eq. (25), that 
is, taking into account the 5-function, Do is not essentially 
self-adjoint. In this case, we must And the self-adjoint ex¬ 
tensions of the operator Do corresponding to different types 
of boundary conditions. Such self-adjoint extensions are based 
in boundary conditions at the origin and conditions at infin¬ 
ity [65-67]. Erom the theory of symmetric operators, it is 
a well-known fact that the symmetric radial operator Do is 
essentially self-adjoint if 

|m + 0 |>l, (28) 

while for 

|m + 0|<l, (29) 


it admits an one-parameter family of self-adjoint extensions 
[33], Dq j-i , where is the self-adjoint extension parame¬ 
ter. According to Ref. [62], the operator (26) admits an one- 
parameter family of self-adjoint extensions. To characterize 
this parameter family of self-adjoint extension, we use the 
approaches proposed by Kay-Studer (KS) [68] and Bulla- 
Gesztesy (BG) [69], being both based on boundary condi¬ 
tions. In short, in the KS approach, the boundary condition 
is a match of the logarithmic derivatives of the zero-energy 
solutions for Eq. (24) and the solutions for the problem Dq 
plus self-adjoint extension. In the BG approach, however, 
the boundary condition is a mathematical limit allowing di¬ 
vergent solutions for the operator Do in Eq. (26) at iso¬ 
lated points, provided they remain square integrable. Then, 
following Ref. [62], the energy spectrum using the KS ap¬ 
proach is found to be 


Po 

+ph 


0 s' 




- 01 ) 


0 s' — \m + <j)\J E(1 — |ot-|-0 |) 


1 

|/n+0| 


(30) 


where po is a finite very small radius. This radius may be 
understood as a kind of physical regularization for the 5- 
function in Eq. (25). Moreover, according to the BG method, 
the energy spectrum is given by 


{E‘mf-M^ = -4 


1 r(l + |m-b0|)\ 
C',r(l-|m-f0|)J 


1 

|m+0| 


+ pI (31) 


r 
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We can note in Eq. (30) that the KS method gives us energy 
levels without any arbitrary parameter that can come from it, 
while the BG method gives an expression that leaves an ar¬ 
bitrary parameter, namely, the self-adjoint extension of the 
parameter However, if we directly compare Eqs. (30) 
and (31), an expression for the self-adjoint extension pa¬ 
rameter is exactly found, i.e.. 


Cm 


2 


-Po 


|m+(/)| 


/0s' - |m-|-0|\ 
\ 0s' -f |m -I- 01 y 


(32) 


The component of the DKP spinor for bound-state solutions 
for the spin projection s* '^ = 1,-1 is given by 


the appropriate boundary condition. In this sense, the self- 
adjoint extension approach was used to determine the bound 
states of vector bosons for the spin projection s^P = 1 , — 1 in 
terms of the physics of the problem, in a very consistent way 
and without any arbitrary parameter. Einally, expressions for 
the bound states energy for vector bosons in the presence of 
the AB potential has been obtained and the conditions in 
which they occur were established. 
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